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We present a theoretical study of a system with competing short-range ferromagnetic attraction and a long- 
range anti-ferromagnetic repulsion, in the presence of a uniform external magnetic field. The interplay between 
these interactions, at sufficiently low temperature, leads to the self-tuning of the magnetization to a value which 
triggers phase coexistence, even in the presence of the external field. The investigation of this phenomenon is 
performed using a Ginzburg-Landau functional in the limit of an infinite number of order parameter components 
(large A' model). The scalar version of the model is expected to describe the phase separation taking place on 
a cell surface when this is immersed in a uniform concentration of chemical stimulant. A phase diagram is 
obtained as function of the external field and the intensity of the long-range repulsion. The time evolution of 
order parameter and of the structure factor in a relaxation process are studied in different regions of the phase 
diagram. 
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I. INTRODUCTION 

Several natural and social processes are governed by com- 
peting interactions. Often the interplay between opposite ac- 
tions produces ordered phases and symmetry breaking events. 
For example, models based on competing interactions are able 
to explain lamellar phases in charged colloids 1 1], pattern for- 
mation in magnetic films, Languimir monolayers and liquid 
crystals ^ and some market behaviors ^ . In living organ- 
isms there is an important example of these processes: the 
spatial orientation of eukaryotic cells 0] called eukaryotic di- 
rectional sensing. Many eukaryotic cells are able to orient 
(polarize) for moving along directions after an external stim- 
ulation. This process is fundamental for important biological 
functions like morphogenesis of organs and tissues, wound 
healing, immune response, social behaviors of some amoe- 
boid cells. The process of orientation takes place on the cell 
membrane where the pattern formation of domains of two dif- 
ferent enzymes determines a symmetry breaking which trig- 
gers the directional sensing jst] . Pattern formation occurs as a 
response to an external stimulation, usually a chemical signal 
activating specific receptors on the cell surface, enhanced by a 
cascade of chemical reactions leading to the cell polarization 
(see Ref. 16i] and references therein). Experimental observa- 
tions 01 suggest that the domain formation is a consequence 
of self-organization of molecular patches. 

Let us briefly summarize the biological mechanism of di- 
rectional sensing. It can be explained in term of the inter- 
play between two enzymes: PTEN (phosphatase and tesin ho- 
molog) and PI3K (phosphatidylinotisol 3-kinase). This inter- 
play is mediated by two lipids: the PIP2 (phosphatidylinoti- 
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sol bisphosphate) and the PIP3 (phosphatidylinotisol trispho- 
sphate). Before stimulation, the cell membrane is populated 
only by the PTEN enzyme with its product PIP2 but, when 
the external chemoattractant is switched on, the enzyme PI3K 
goes from cytoplasm to the cell membrane and binds to re- 
ceptors. Then, the interplay between the two enzymes takes 
place: the PI3K catalyzes PIP2 in PIP3 and PTEN catalyzes 
PIP3 in PIP2. The enzymes can bind to the respective lipid 
products which diffuse over the membrane. Enzymes can un- 
bind from membrane and quickly diffuse in cytoplasm binding 
again in another place of the membrane. Catalysis and lipides 
diffusion mediate an effective short-range attraction between 
enzymes of the same type. The quick diffusion of enzymes 
in the cytoplasm mediates a long range interaction. The com- 
bination of these actions produces the phase separation of a 
PI3K rich zone and a PTEN rich zone. 

The natural framework to treat the process, from a physi- 
cal point of view, is the statistical physics of phase separation, 
which is useful to understand and describe in synthetic way 
the behavior of many complex systems that give rise to orga- 
nization and pattern formation f8^,xl- The spatial organization 
phenomena described above have the characteristic of self- 
organized phase separation processes, where the cell state, 
driven by an external field, decays into the coexistence of two 
chemical phases, spatially localized in different regions. It 
was shown recently, by Monte Carlo simulation in a lattice 
gas model [10], that the phenomenology of directional sens- 
ing can be obtained by using an effective free energy. A simi- 
lar point of view was used in the recent papers of Gamba et al. 
Ill L 121. The remarkable physical characteristic of the process 
is that the orientation is possible for a wide range of external 
chemical attractant flS*]. Namely, the phase coexistence and 
separation is possible for different values of an external field. 
In the lattice gas model Holl . a long range repulsion, derived 
from the interaction with a finite cytosolic reservoir of total 
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enzymes and the interaction with an external field, modelling 
the action of the chemical attractant, give rise to coexistence 
for an interval of values of the external field. A short-range 
attraction between enzymes, derived from their catalytic ac- 
tions on lipids and diffusion, gives rise to a coarsening pro- 
cess which produces phase separation. A similar mechanism 
operates in some econophysics models based on two major 
conflicting interactions in economy, the tendency of a trader 
to follow the actions of his neighbors and the tendency to fol- 
low the actions of the minority. In the language of magnetic 
systems, which we will use throughout this paper, the first 
tendency can be modelled by a short-range ferromagnetic in- 
teraction, while the second one by an antiferromagnetic long 
range interaction. 

Inspired by the mechanism of eukaryotic directional sens- 
ing and motivated by its wide applicability, here we present 
the analytical treatment of a system with competing short- 
range attraction and long-range repulsion, under the action of 
a uniform external field. This is done in the framework of the 
time dependent Ginzburg-Landau (TDGL) theory for the evo- 
lution of the order parameter. The analytical tractability of the 
dynamics is achieved considering a vector order parameter in 
the limit of an infinite number of components (large limit) 
ifTill . The equations of motion are derived in the scheme of 
non conserved order parameter, corresponding to the absence 
of local concentration of enzymes. 

The large limit is a powerful method to obtain analyti- 
cal results, nevertheless it is necessary to keep in mind some 
important differences with the nonlinear models usually em- 
ployed in the description of phase ordering when the order pa- 
rameter is a scalar, such as kinetic Ising models or the TDGL 
with (|)^ interaction 1 15]. The basic difference is in the mecha- 
nism of equilibration after a symmetric quench below the crit- 
ical point. In the nonlinear models the system responds to the 
dynamical instability by the formation and growth, through 
coarsening, of domains of the ordered phases. This leads to 
the development of a bimodal probability distribution for the 
local magnetization, which eventually ought to evolve into the 
symmetric mixture of the two possible broken symmetry or- 
dered phases llT6ll . This we call an ordering process. In the 
large limit, instead, the development of a bimodal distribu- 
tion, and therefore ordering, is not possible, since, as it will be 
clear below, the system is effectively linearized and the statis- 
tics are Gaussian. Then, the response to the dynamical in- 
stability takes place through the development of macroscopic 
fluctuations in the most unstable Fourier component of the or- 
der parameter, through a process which is formally identical 
to the one leading to the formation of the condensate in the 
low temperature phase of the ideal Bose gas. We refer to this 
equilibration mechanism as condensation of fluctuations. The 
differences between the two equilibration processes have been 
investigated in detail in Ref. [17, 18]. 

The remarkable feature of the large limit, and the reason 
for its wide use, is that despite the considerable difference in 
the physical processes of equilibration, the phenomenology of 
the observables of interest, such as correlation functions and 
response functions, is the same as in the nonlinear models, 
apart for obvious quantitative discrepancies, like the values of 



exponents or the shape of scaling functions. The typical ex- 
ample is that of the equal time structure factor, which displays 
dynamical scaling and the growth of the Bragg peak in the 
large A^ limit |19], exactly as in the nonlinear models. It is, 
then, matter of interpretation in one case to read the growth 
of the Bragg peak as revealing domain coarsening and, in the 
other, condensation of fluctuations. There is, by now, a vast 
body of literature documenting the robustness of the large A^ 
limit in reproducing the phenomenology of phase ordering in 
a large variety of models, warranting to overlook the distinc- 
tion between ordering and condensation, as we shall do in this 
paper, when one is interested in the main qualitative features 
of the process. 

The paper is organized as follows. In section II we carry out 
the large A^ limit on the TDGL model for a vector order pa- 
rameter, deriving the basic equations. In section III we study 
the equilibrium properties of the system, obtaining the phase 
diagram in the temperature and external field plane, parame- 
terized by the strength of the long-range repulsion. This de- 
limits the region of parameters where condensation or, equiva- 
lently, phase coexistence and separation are possible. Section 
IV is devoted to the study the time behavior of the average 
value of the order parameter (magnetization) and of the cor- 
relation function. We recall that in the scalar case the mag- 
netization represents the concentration difference between the 
two species of enzymes. Concluding remarks are presented in 
section V. 



II. THE MODEL 

The system is modelled by a free energy functional of the 
form 
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where (|)=((|)i,...,(|)Af)isanA^ component vector order param- 
eter We shall take r < 0, ^ > 0, i/ ~ O (A^'/^), V is the volume 
of the system and A- > is an antiferromagnetic coupling. We 
shall consider the static and dynamic properties of the model. 

As is well known from the theory of critical phenomena, the 
introduction of an A^ component order parameter is a very con- 
venient technical device to generate controlled and systematic 
correction about mean field behavior, using 1 /N as an expan- 
sion parameter [201 . We shall limit the treatment to the lowest 
(mean field) order by taking the large A^ limit (A^ — > 



A. Equation of motion 

In the framework of TDGL model for the dynamics 
3(j)(x,f) 5yf[(j)] 



dt 



{x,t)+f[{x,t) 



(2) 



3 



where r[{x,t) is the white noise at temperature T, with zero 
average and correlator 

the equation of motion of the order parameter is given by 



dt 



X, f ) + r f ) + — ((|) ■ (|))(|) (x, f ) - H(x) 



N 



— / dx(^{x,t)+r\{x,t). 
V Jv 
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Since the external field breaks the rotational symmetry in the 
order parameter space, it is convenient to introduce the longi- 
tudinal and transverse components with respect H{x) 

4 = ^\\+^± (4) 
and then to split the longitudinal component into the sum 

$11 (x,f) =M(x,f)+ V|/(x,f) (5) 

where M(x,f) = ($||(x,f)) is the magnetization, while the av- 
erage longitudinal fluctuations vanish (\|/(x,f)} = by con- 
struction. The angular brackets denote the average over both 
the initial condition and the thermal noise. In the following 
we shall take a reference frame with the 1-axis along the lon- 
gitudinal direction. 

Assuming, next, M ^ o{N^I^), xj/ ~ 0(1) and comparing 
terms of the same order of magnitude in A^, to leading order 
we get the pair of equations 
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\—V^ + r + gm^+gS\m+(h / dxm) 

V Jv 
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and 
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where m = m{x,t) and h = h{x,t) are the following rescaled 
quantities 

m{x,t)^M{x,t)/N^I^, h{x,t)=H{x,t)/N^/^. (8) 

In the large limit the quantity S{x,t) is given by the self- 
averaged fluctuations 
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of the generic transverse component (])«. Furthermore, since 
in Eq. Q the components of (j)^ are effectively decoupled, 
from now on we shall refer to the equation for the generic 
component omitting the a subscript. 

Taking a space and time independent external field h and 
space translation invariant initial conditions, we can assume 
space translation invariance to hold at all times. Hence, 
Fourier transforming with respect to space, and introducing 
the equal-time transverse structure factor 



we obtain the closed set of equations 

dm{t) 



dt 



dC±(k,t) 
dt 



(o(p,t)m{t) + h 



= -2(o{k,t)C^{k,t)+2T 



(11) 
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with m(f) = m(k = 0,?) and &(k,t) defined by 

(a{k,t)=k^+Xdk,o + r + g{m^ + S) 
and the noise correlator in Fourier space given by 



With periodic boundary conditions the wave vector runs over 
k = ^n, where n is a vector with integer components and 

L'' = y . Furthermore, sums over k like the one in Eq. (fT3l i are 
cutoff to the upper value k,„ax = A, where A^^ is related to a 
characteristic microscopic length, for instance the lattice spac- 
ing of an underlying lattice. Finally, the longitudinal fluctua- 
tions \\t have been dropped since do not give any contribution 
to leading order. 



III. STATIC PROPERTIES AND PHASE DIAGRAM 

If equilibrium is reached, all quantities become time inde- 
pendent. Rewriting Eq. (fT4l i as 

f X + jj, for ^ = ^, 
«(^) = i^2/^^ ,^0 (15) 



with 



^l = r + g{m^ + S) (16) 



and putting to zero the time derivatives, from Eqs. (fTTT i. (fT2l i 
and (fTJl i we obtain the set of equations 



(k + ij)m = h 



(a{k)C^{k)^T 



(17) 



(18) 



(19) 



From Eqs. ([TOl) and ^ foUows C± (k) > and (x){k) > 0, re- 
spectively. The latter inequality, because of Eq. ST5[ . requires 
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where ^„„„ ~ 1/L is the minimum allowed value of ^ ^ 0. 
Therefore, for a given A- and for V sufficiently large. A- +/j > 
and Eq. ( fTTI i can be rewritten as 



(21) 



In the same way, from Eq. ( fTSI l we can write 



for 
for 



k = Q 
k > k„ 



(22) 



where C±(k,„i„) diverges as /j approaches /j„„>,. Notice that /v 
can be identified with the inverse square transverse correlation 
length 

In order to obtain the full solution, we must now determine 
how jj depends on the parameters of the problem {T,h,V,X). 
In principle, this can be done by inserting the above results 
into Eq. ( fT6b and solving the basic self-consistency equation 



VX- 



Ily ^ 

V P k'-^ 

kytO 



A* 



(23) 



However, for our purposes general considerations are suffi- 
cient, without actually solving the above equation. For V suf- 
ficiently large the second term in the right hand side can be 
neglected and the equation can be rewritten in the form 



X- 



V fJ-lLl,, 



V -4- 



1 



k^- 



(24) 



where the A;„„„ term has been extracted from underneath the 
sum. Letting /j to vary from /j„„>, to qq, the left hand side 
is a monotonously increasing function of /j, while the right 
hand side diverges at /j„„„ and decreases monotonously with 
increasing /j. Therefore, for any finite V, there exists a solu- 
tion ifiy) > limin- Looking at Eqs. (l2Tl i and (|22] |. this means 
that the system behaves paramagnetic all over the (T, /z) plane, 
with a finite structure factor The difference with respect to 
what one would have in the purely short-range model, due to 
A, ^ 0, is revealed by the anomaly (l22b in the structure factor 
at ^ = and by the reduction of the magnetization in Eq. ( |2TI) . 
Rewriting the latter as /j* (V )m = heff with 



h — Xm, 



(25) 



we see that the reduction of the magnetization comes about 
through a feedback mechanism, whereby the external field h, 
via the long-range interaction, is substituted by h^ff. 

Let us now see what happens in the infinite volume limit. 
From Eq. ( |20| ) follows /j„„„ 0^ and there are two possibili- 
ties 



lim fj*{V) 



IJ*>0 
IJ*=0. 



(26) 



In the first case, the second term in the right hand side of 
Eq. ( l24b can be neglected yielding 



(27) 




FIG. 1 : (Color online) Phase diagram showing the X dependence of 
the region within which condensation of the transverse fluctuation 
takes place. The presence of a condensate, in the blue regions, cor- 
responds to phase separation and coexistence in the corresponding 
scalar model. In this and in all the other figures g = —r = \, A = 2ti 
and rf = 3, yielding 7c (0) = 71. 



where 



B(u) = hm - y 



(28) 



is a monotonously decreasing function of /j with the maximum 
value 



B{0) 



d-2 ' 



for 



d>2 



divergent for d <2 



(29) 



with Kd = [2^"%''/2r(<i/2)]"' being the <i-dimensional solid 
angle [21J and Y{d/2) the Euler gamma function. Hence, 
Eq. dZTb admits a positive solution if T is greater than the h- 
dependent critical temperature 



Tc{h) 



+ 8hVX^ 
gBiO) 



(30) 



which vanishes for d <2 and for any h, while it is finite for 
d > 2 reaching the maximum value Tc{0) = —r/gB{0) for 
h = and decreasing to zero when h reaches the limit values 
±hc (see Fig. [T]i with 



hc = X{~r/g)'/\ 



(31) 



Conversely, if T and h are such that T < Tc{h), Eq. (|27] | 
cannot be satisfied. This means that the second of the two pos- 
sibilities in Eq. ( |26] i applies, requiring to keep also the second 
term in the right hand side of Eq. (|24] |. which we rewrite as 



A* 



—C±ik„ 



Since this is satisfied for ^ = 0, we get 



C±{kmi„) =V 



Tc{h)-T 



Tc{h) 



(32) 



(33) 
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showing that C±{kmi„), for T < Tc{h), diverges like the vol- 
ume in order to give a finite contribution in the right hand side 
ofEq. (El. 

Summarizing, in the infinite volume limit 



^ r h/{X + p*), for T>Tc{h) 
\h/X, for T <Tc{h) 



(34) 



where /j* > 0, and 



r m = / ^/(^'+^*)' T>Tc{h) 

' \T/k^ + 'M^{T)h{k-Q+), for T<Tc{h) 

(35) 

showing that there is condensation of transverse fluctuation at 
A: = 0+ for T < Tc{h), with the size of the condensate given 
by 



Tc{h) 



Tc{h) 



(36) 



In order to understand this result, it should be recalled that 
in the purely short-range large model the phase transition 
occurs only on the h = Q axis, where for T < 7c (0) the con- 
densation of fluctuations takes place (v^ at k = 0. Conden- 
sation of fluctuations means that C± (k = 0) becomes macro- 
scopic in order to equilibrate the system below Tc{0) with- 
out breaking the symmetry, through a mechanism very simi- 
lar to that of the Bose-Einstein condensation, as mentioned in 
the Introduction. No other mechanism is available, since the 
large limit renders the system effectively Gaussian iflTlfTsll . 
However, condensation of fluctuations produces the onset of 
a Bragg peak at A: = and, therefore, a phenomenology of the 
structure factor which is indistinguishable from that due to the 
occurrence of phase separation in the nonlinear models ifTsll . 
When li y^Q the symmetry is broken and equilibrium can be 
established at any temperature through the development of a 
non vanishing magnetization, without any condensation. 

In the system with the long-range coupling of antiferromag- 
netic type everything remains the same along the h = axis, 
since the symmetry is unbroken, the magnetization is zero and 
the only effect of the X term is to shift the Bragg peak from 
^ = to ^ = 0+. The novelty appears outside of the h — Q 
axis, where the explicitly broken symmetry induces the de- 
velopment of a non vanishing magnetization which, however, 
through the feedback mechanism driven by the antiferromag- 
netic interaction produces the effective reduction ( |25] ) of the 
external field. So, if for a given X, the values of T and h man- 
age to make h^ff = 0, at that point the value of the magne- 
tization gets stabilized to the value (|34] | and the only way to 
equilibrate the system is through the condensation of fluctua- 
tions. The result is the phase diagram of Fig. [T] showing the 
expansion of the phase coexistence region outside the li = 
axis. The constant X curves delimit the regions on the {T,h) 
plane within which the system self-tunes the final magnetiza- 
tion to the value h/X such that heff = triggering, therefore, 
the condensation of the k — Q^ fluctuations. 
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FIG. 2: (Color online) Time evolution of the magnetization for dif- 
ferent values of h,X= I and T = 0.25 Tc{h). 



IV. DYNAMICAL PROPERTIES 

In order to investigate the dynamics, we have solved nu- 
merically the coupled equations (fTTI) . (fTSl i and ( fTSI ) in a dis- 
cretized three-dimensional Fourier space using fourth-order 
Runge-Kutta method with adaptive step size f2Z]. We have 
used a mesh of linear size L — 1000, taking A- = 1 and the 
initial conditions m(0) ~ —1, C^(fe,0) = 0. 

Let us first consider (T, /z) in the region of phase separation, 
with T < Tc{h). Fig.|2]and Fig.[3]illustrate quite well the con- 
siderations made at the end of the previous section. The first 
one displays the evolution of the magnetization for different 
values of h. After a fast transient there is saturation to the 
equilibrium value 



h 



meq = ^ 



(37) 



taking place from above or from below for h positive or nega- 
tive, respectively. Apart from the details of the transient, this 
behavior of the magnetization agrees with the prediction and 
results of Ref. 1 10] for the scalar model. 

As explained above, when the magnetization gets stabilized 
at the value ( [37l ) by heff = 0, the growth of the Bragg peak 
is inevitable in order to equilibrate the system. This is illus- 
trated in Fig. [3] In the late time regime the structure factor is 
expected to obey a scaling form of the type [.15il 



C^{k,t)^L''(t)F{kL{t)), 



(38) 



where F{kL(t)) is a scaling function and L(f) ^ f'/^ with 
z = 2, as appropriate for phase ordering processes without 
conservation of the order parameter, is a time dependent char- 
acteristic length |23]. The same growth law f'/^ was found in 
Ref 111 01 for the mean cluster size. The inset of Fig.[3]shows 
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FIG. 3: (Color online) Time evolution of the transverse structure fac- 
tor C^(k,t) for h = 0.5, X = 1 and T = 0.25Tc{h). Inset: time evo- 
lution of the peak height C±_ {k„,i„ ,l) 



FIG. 4: (Color online) Behavior of magnetization outside the coex- 
istence region, with h = 0.5, A. = 1 and T = 2Tc(h). The dotted line 
represents the equilibrium value of magnetization for T < Tc{h). In- 
set: time evolution of C±{k,„i„,t). 



that the height of the peak follows quite well the power law 

C(^™„,0~f^/^. (39) 

Actually, for times of order 10^ — 10'* there appears a devi- 
ation from the power law ( [39] l. This is a finite size effect, 
unavoidable in the numerical computation and not to be con- 
fused with the equilibration of the magnetization, which is in- 
dependent of the size of the system (notice the huge differ- 
ence in the time scales). Therefore, considering the infinite 
system, we have the interesting instance of two observables 
in the same system, one of which (the magnetization) equi- 
librates rather quickly, while the other (the structure factor) 
does not reach equilibrium in any time scale. 

Conversely, if T > Tc{h) there is a solution of the self- 
consistency relation without growth of the condensate and 
both the magnetization and the structure factor reach equilib- 
rium in the same time scale. This is illustrated in Fig. |4] dis- 
playing the saturation of the magnetization to the limit value 
nieij — with /J* > 0, in agreement with Eq. ( l34l i. while the 
inset shows the saturation of the peak height to the equilibrium 
value 

T/{kl„+iL,*) (40) 

in agreement with Eq. ( l35T l. 

The same qualitative behavior, that we have here illustrated 
in the d ~3 case, is expected for any d > 2. For d — 2, due 
to the divergence in the denominator of Eq. dSOl l. the critical 
temperature vanishes squeezing the coexistence region of Fig. 
[T]onto the vertical axis at T = 0. Therefore, in the d — 2 case 
condensation of fluctuations can be obtained only for T = 
and < he- 



V. SUMMARY 

In this paper we have studied the static and dynamic prop- 
erties of a system described by a free energy functional with 
a short-range ferromagnetic interaction and a long-range anti- 
ferromagnetic interaction, in presence of an external uniform 
magnetic field. The analysis has been carried out in the large 

limit. The scalar counterpart of this model is the lattice-gas 
Hamiltonian 1 10] used to model the phenomenology of phase 
separation occurring in the inner part of cell surface during 
directional sensing. We have focussed on the phase ordering 
process taking place below the critical temperature, even in 
the in presence of the external magnetic field. 

In particular, through the equations of motion for the mag- 
netization and the transverse structure factor, we have high- 
lighted how the competing interactions induce the self-tuning 
of the magnetization within the phase coexistence region. We 
have derived the phase diagram, which depends on the mag- 
netic field and the strength of the antiferromagnetic coupling, 
showing that phase-separation is possible for a range of val- 
ues of the external field. Taking the large limit it has been 
possible to derive analytically the dependence of the critical 
temperature on the magnetic field and on antiferromagnetic 
coupling X. The phase diagram of Fig. [T]depicts in the {T,h) 
plane the phase coexistence regions for different values of X. 
The equilibrium value of magnetization nieq —h/Xis in agree- 
ment with the prediction and with the results obtained in the 
Monte Carlo simulation for a lattice gas system (Ref. [10]). 
The dynamics shows that the antiferromagnetic coupling com- 
bines with the magnetization to generate the effective mag- 
netic heff field, eventually vanishing within the coexistence 
region. While the magnetization equilibrates very quickly, the 
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structure factor does not equilibrate on any time scale. 

For T < Tc{h) the relaxation process is characterized by 
the growth of a condensate in transverse structure factor at the 
most unstable wave vector k = kmi„. The onset of condensa- 
tion signals the occurrence of a phase separation and corre- 
sponds to domain coarsening in the scalar case. The late time 
behavior of the structure factor is characterized by dynamical 
scaling and power law growth of the peak, with the exponent 
d/2 characteristic of non conserved order parameter. 

In conclusion, we have analyzed a model where the com- 
petition between the short-range ferromagnetic and the long- 
range antiferromagnetic interaction between two species leads 
to the phase separation for a wide range of external field and 
temperature. The occurrence of phase separation is a crucial 



intermediate step allowing for the amplification of external 
field gradients, leading to directional sensing as illustrated in 
Ref. U^. 

The general property of the free energy functional ([T]l, giv- 
ing rise to phase coexistence through self-tuning, can be very 
useful in other contexts characterized by the balance between 
short range attraction, long-range repulsion and an overall ex- 
ternal action. 
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